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LQR  for  in  situ  discrete  structural  damage  growth  retardation 

Chimpalthradi  R.  Ashokkumar1  and  George  WP  York2 
US  Air  Force  Academy,  2354  Fairchild  Dr,  Ste  2F6,  Colorado  Springs.  CO  80840-6236. 

ABSTRACT 

One  of  the  objectives  in  structural  health  monitoring  is  damage  localization  and  assessment. 
During  service  intervals,  damages  are  localized;  however,  assessment  (the  sizes  of  damages)  is 
generally  difficult  to  determine  and  it  requires  the  in  situ  data.  Transient  and/or  steady  state 
vibration  database  (real  time  data  or  offline  data  in  frames)  such  as  displacements,  velocities 
and/or  accelerations  is  a  useful  resource  for  this  purpose.  Yet  growth  assessment  remains  to  be  a 
challenging  problem.  Growth  retardation  is  possible  when  the  tensile  and  compressive  action  at 
the  damaged  spring  is  minimized.  Hence  the  problem  is  posed  as  follows.  When  unknown 
damage  is  present  in  one  of  the  linear  springs  of  a  discrete  structure,  in  this  paper,  in  situ 
localization  is  pursued  and  then  its  growth  is  mitigated  by  applying  the  linear  quadratic  regulator 
(LQR)  design.  Although  the  technique  applies  to  multiple  damages  with  an  unknown  structural 
load,  in  this  paper  only  single  damage  and  transient  vibrations  arc  considered.  Another  important 
requirement  for  damage  growth  retardation  is  to  keep  the  real  parts  of  the  eigenvalue 
distributions  of  the  damaged  structure  farther  at  the  left  half  plane  of  the  complex  plane 
compared  to  the  pristine  structure.  This  requirement  depicts  robust  alpha  stability  problem.  A 
discrete  structure  is  considered  for  illustrating  several  of  these  features  considered  in  this  paper. 

Key  words :  In  situ  damage  detection,  growth  retardation,  vibration  data,  robust  alpha  stability. 
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I,  Introduction 


A  wide  range  of  structural  health  monitoring  (SHM)  problems  focuses  on  damage  detection 
[1-3],  localization  [4-7],  assessment  [8,9],  and  mitigation  [10,11]  under  various  loading 
conditions.  A  structure  with  damage  is  often  detected  during  service  intervals  in  offline  although 
the  damage  occurs  in  real-time.  Thus  from  the  time  the  damage  is  occurred  to  the  time  the 
damage  is  detected,  a  structural  failure  with  damage  growth  can  happen.  In  this  context,  in  situ 
damage  localization  and  damage  growth  retardation  becomes  an  important  problem.  It  helps  to 
extend  the  remaining  service  life  of  the  structure.  For  damage  growth  retardation,  note  that 
model  prediction  by  damage  assessment  is  required  to  determine  an  actuator  load  and  mitigate 
damage  growth  online  at  least  until  a  repairing  or  a  replacement  is  performed.  Hence,  damage 
localization  and  damage  growth  retardation  continue  to  be  one  of  the  important  SHM  problems. 
Consider  a  discrete  dynamic  structure  made  of  several  springs,  masses,  and  damper  components 
operating  with  a  vibration  (transient  or  steady-state)  database.  These  systems,  for  instance,  may 
represent  the  flexible  multi  body  belt  drive  system  modeled  using  a  technique  such  as  the 
Lagrangian  approach  [12].  In  large  scale  systems,  the  structure  may  consist  of  several  springs 
and  if  one  of  the  springs  is  damaged,  then  the  whole  system  may  lead  to  a  failure.  It  is  also 
important  to  check  periodically  the  health  of  each  spring.  In  these  scenarios,  assuming  that  one 
of  the  linear  springs  of  the  structure  is  required  to  be  monitored  for  possible  damage,  new 
procedures,  based  on  the  fact  that  the  reduced  stiffness  at  the  damaged  spring  exhibit  more 
displacement  activity,  are  captured  using  hat  functions  defined  for  that  spring  [10].  In  this  paper, 
in  situ  damage  localization  is  pursued  by  applying  this  technique  and  after  the  damage  is 
localized,  damage  growth  retardation  is  pursued  using  the  LQR  formulations  [13].  Vibration 
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database  is  utilized  for  damage  localization.  In  retardation,  a  fact  based  on  the  robust  alpha 
stability  criterion  to  restrict  the  eigenvalue  distributions  of  the  damaged  structure,  whether 
controlled  or  uncontrolled,  is  applied  to  limit  the  damage  growth  beyond  which  the  structure  may 
vibrate  worse  than  the  pristine  structure.  Here  alpha  refers  to  the  real  part  of  the  eigenvalue  of  the 
pristine  structure  that  is  nearest  to  the  imaginary  axis. 

Model  updating  for  damage  growth  retardation  is  possible  when  stiffness  loss  at  the  spring  is 
known  as  an  indication  of  damage.  Stiffness  loss  due  to  damage  is  extensively  analyzed  in 
composite  structures  [14-18].  If  finite  elements  are  used,  stiffness  loss  introduces  perturbations  to 
the  finite  element  stiffness  matrix  distributions  across  the  continuous  structure.  In  discrete 
structures,  depending  upon  the  damaged  spring,  the  stiffness  matrix  will  exhibit  a  distribution 
matrix  at  which  the  stiffness  of  the  spring  will  get  deteriorated.  The  vibration  database  is  utilized 
to  determine  the  distribution  matrix  [10].  Hence  the  model  is  updated  with  the  damage  size.  The 
hat  functions  considered  in  this  paper  to  determine  spring  stiffness  loss  in  the  stiffness  matrix 
also  applies  to  continuous  structures  if  vibration  based  SHM  techniques  pursued  in  the  literature 
[19-21]  suggest  various  uncertain  finite  element  stiffness  loss  distributions  across  the  continuous 
structure.  Yet,  a  disadvantage  with  the  present  method  in  [10],  even  for  discrete  structures,  is  that 
displacements  at  all  locations  of  the  masses  are  assumed  available.  In  state  feedback  format,  the 
technique  may  be  pursued  further  by  estimating  all  the  states  variables  using  an  observer.  One  of 
the  merits  of  this  approach  is  that  the  damage  localization  technique  is  numerically  attractive. 
However,  the  pristine  stiffness  data  of  the  springs  are  assumed  to  be  known  perfectly.  In  practice, 
when  the  structure  becomes  old,  the  stiffness  of  the  linear  springs  could  have  been  deteriorated. 
In  this  case,  multiple  damage  localization  becomes  important.  Hence  the  approach  presented  in 
this  paper  does  accommodate  these  scenarios  of  the  SHM  [10].  Although  the  procedure 
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presented  in  this  paper  is  extendible  to  multiple  springs,  presently,  damage  localization  and 
retardation  is  pursued  when  one  spring  at  a  time  is  considered  to  be  damaged. 

After  the  damaged  spring  is  identified,  the  second  objective  of  this  paper  is  to  retard  the 
damage  growth  so  that  the  amount  of  stiffness  loss  in  the  damaged  spring  is  minimized  during  its 
operation.  The  hat  functions  derived  previously  for  in  situ  damage  detection  are  further  utilized 
to  minimize  the  tension  and  compression  activity  (can  be  defined  as  transient  fatigue  cycles 
similar  to  the  structural  responses  resulting  from  a  fatigue  load).  Although  LQR  is  suitable  to 
handle  hat  functions  associated  with  the  tensile  and  compressive  forces  (contributing  to  potential 
energy)  and  their  rates  (transient  fatigue  rate  cycles  contributing  to  the  kinetic  energy)  so  that  the 
activity  at  the  damaged  spring  is  minimized,  the  question  arises  on  the  linear  model  to  be  used  to 
design  LQR.  As  long  as  the  unknown  eigenvalues  of  the  closed  loop  system  matrix  for  the 
damaged  model  are  farther  to  the  left  half  plane  of  the  complex  plane  compared  to  the  pristine 
model’s  eigenvalues,  it  is  known  that  transients  of  the  controlled  structure  with  damage  decay  at 
faster  rate  than  the  uncontrolled  pristine  structure.  Therefore,  it  is  the  real  eigenvalue  of  the 
damaged  structure  that  must  be  kept  farther  from  the  eigenvalue  of  the  pristine  structure  which  is 
nearest  to  the  imaginary  axis.  Again  damped  frequency  (imaginary  part  of  the  uncertain 
eigenvalue  hopefully  residing  in  open  left  half  plane  of  the  complex  plane)  of  oscillations 
contributes  to  damage  growth  retardation.  Presently  only  real  part  restrictions  are  imposed.  This 
investigation  leads  to  the  determination  of  potential  damage  sizes  that  will  make  the  eigenvalues 
of  the  damaged  structure  to  cross  alpha  limit  which  is  the  eigenvalue  of  the  pristine  structure 
nearest  to  the  imaginary  axis.  In  this  paper,  these  damage  sizes  are  presented  using  robust 
stability  analysis  [22],  In  real  time,  techniques  using  EKF  are  reported  in  the  literature  [23-25]  to 
estimate  these  parameters.  When  the  model  is  updated  with  the  potential  stiffness  loss  parameters 
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[26,27],  actuator  loads  to  mitigate  damage  using  complete  or  partial  pole  placement  or  adaptive 
control  are  presented  [28,29].  In  this  paper,  LQR  is  performed  to  minimize  potential  and/or 
kinetic  energy  type  of  functions  when  hat  functions  are  utilized.  Then,  using  alpha  stability 
criterion,  admissible  damages  for  damage  growth  retardation  are  presented. 

The  paper  is  organized  as  follows.  In  Section  2,  in  situ  structural  damage  localization 
problem  is  considered.  After  the  damage  is  localized,  in  Section  3,  damage  growth  retardation 
using  LQR  design  is  pursued  by  using  the  pristine  model.  Then  admissible  damages  retarding  the 
damage  growth  are  presented  in  Section  4.  Simulations  and  conclusions  are  presented  in  Sections 
5  and  6,  respectively. 

2.  Discrete  In  Situ  Structural  Damage  Localization 

In  this  section,  a  problem  formulation  for  in  situ  damage  localization  using  a  vibration 
database  available  in  real  time  as  a  raw  data  or  as  a  stored  database  in  frames  is  presented. 
Consider  a  second  order  spring-mass-damper  system  representing  the  mechanical  systems  such 
as  the  one  considered  in  [10,1 1,28]  and  shown  in  Figure  1  as  follows. 


g— m— 

e, 


- - 

— *<h 

1 

! 

J 

_  *3  _ 

m, 


/Wl 

— Q—  w- 

e , 


— m— 
— D— 

e~ 


tth 


*ch 

kl  -  Ak, 

—m — | 


Figure  1  A  Discrete  Structure  Representing  a  Mechanical  System  for  Damage  Localization  in 

Either  Spring  2  or  Spring  4 

M$(0  +  C0q(0  +  K0q(/)  =  0,  q(Q)*0  (1) 

where  q(/)e/?"  is  a  vector  containing  n  damage-free  structural  degrees  of  freedom;  4(0  and 
q{0  respectively  represent  the  velocity  and  acceleration  of  q{/) .  M,  K0  and  Co  are  mass,  stiffness 
and  damping  matrices  that  are  symmetric  with  order  n  .  Consider  a  linear  spring  at  which  one 
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wants  to  detect  the  damage.  If  a  single  damage  is  present  at  the  /h  spring,  uncertainties  Ak ;  are 
introduced  in  the  stiffness  and  damping  matrices  as  below, 

K  =  K  o  ~  Ak  f  K  j  . 

C  =  C0  -  Ae,C  j . 

K '.j  and  cyare  the  known  distribution  matrices  that  can  be  determined  using  the  stiffness  and 

damping  matrices  K  and  C.  The  second  order  system  in  Eq.  (1)  is  accordingly  modified  as: 

Mij(O  +  Cq(O  +  Kq(/)  =  0,  q(0)*0  (2) 

Note  that  the  states  with  the  damage  are  denoted  by  q(Y).  Damage  detection  (or  localization)  at 
springy  (to  develop  the  technique  assume  the  damper  is  healthy)  using  the  vibration  database  is 
formulated  as  follows.  Assume  q(/)  are  measurable  or  estimated  using  an  observer.  Due  to  the 
stiffness  loss  at  springy,  it  is  expected  that  the  displacements  of  the  masses  between  the  springy 
are  more  for  damaged  system  in  Eq.  (2)  than  the  pristine  system  in  Eq.  (1).  This  behavior  of  the 
damaged  spring  can  be  studied  best  by  using  the  isolated  hat  functions  (IHFs)  as  depicted  in  Fig. 
2.  Let  kj  be  the  springy  between  masses  ma  and  mb.  Consistent  with  the  Figure  1,  the  positive 
directions  of  the  displacements  of  masses  ma  and  mb  are  accordingly  inferred  by  qa(t)  and  qb(t), 
respectively.  IHFs  simply  present  masses  ma  and  mb  in  tension  (or  compression).  For  instance  the 
IHF  in  Fig.  2a  represents  the  damaged  spring  in  compression  with  positive  displacement 
direction  for  qa(t)  and  negative  for  qb(t).  In  the  following  sections,  a  procedure  based  on  these 
IHFs  to  localize  damage  is  presented. 

In  order  to  investigate  displacements  during  compressive  and  tensile  loads  at  springy,  define  the 

modal  vector  for  masses  ma  and  mb  as  [x„,]=7-|p  ■■■  -t  m  o]  (prime  indicates  transpose 

l  „ 

of  the  vector).  Here,  i  v  represents  the  length  of  the  vector  [xm].  I  and  m  are  the  displacement 
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units  for  the  masses  tna  and  mb,  respectively.  The  orthonormal  vector  will  be 
[xm  1  =  •“  m  (  •••  o]  •  This  modal  vector  depicts  the  masses  ma  and  mb  in  phase.  Note  that 


the  modal  vector  |xm]  depicts  the  tensile  load  using  the  hat  function  at  Figure  2b.  The 
compressive  load  is  inferred  as  follows.  Define  an  isolated  modal  matrix  T  as. 


T 


where  £,m  >  0  and  fv 


4?  +  m2  . 


The  displacements  qa(t)  and  qb(t)  are  related  by  the  isolated  modal  matrix  T  as  below, 
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(3) 


b) 


V  ftp* 

Sri 

t  . 

■ — -AM — i 

mb 

m,  — AM —  })] 


Figure  2  IHFs  for  a)  compressive  and  b)  tensile  load  tests  in  springy 
Clearly,  when  ca  is  negative,  the  springy  takes  compressive  load.  Damage  detection  criterion  at 
the  springy  using  the  vibration  database  qjyf)  and  qb(t),  IHFs  and  isolated  modal  matrix  T  is 
presented.  Initial  condition  response  is  employed  to  establish  a  damage  detection  criterion.  When 
the  vibration  database  q(0  due  to  an  initial  condition  is  given,  the  displacements  qa(t)  and  qb(i) 
connected  to  the  springy  through  masses  ma  and  mb,  respectively,  are  known  to  be  measured  or 
estimated  using  an  observer.  Through  the  isolated  modal  matrix  T,  the  constraint  in  Eq.  (3)  is 
established.  Note  that  the  vectors, 
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v  =  —  [-  t  m\  and  w  =  —  \m  t\ , 


(4) 


are  orthonormal  as  in  the  case  of  modal  vectors  for  an  undamped,  uncontrolled  free  vibration 
problem.  For  example,  consider  the  case  n  =  3 .  As  in  the  case  of  a  modal  matrix,  the  vectors  v 
and  w  form  an  orthonormal  basis  for  R3  as  given  below: 
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For  t  >  m ,  tension  and  compression  with  displacement  qiit)  are  considered  for  analysis.  For 
m  >  i ,  tension  and  compression  with  displacement  qi(t)  are  considered  tor  analysis.  In  this  case, 
a  linear  spring  between  the  masses  m2  and  m 3  is  chosen  for  damage  detection.  Similarly,  for  a 
spring  with  displacements  qx(()  and  q2(t)>  following  modal  transformation  can  be  considered  for 
damage  detection. 

(6) 

These  transformation  matrices  can  be  inferred  for  a  generic  displacement  vector  q(/)  with  n 
components.  In  order  to  detect  damage  in  the  spring  j,  vibration  database  for  pristine  and 
damaged  structures  are  assumed  given.  This  is  a  valid  assumption  because  prior  to  damage, 
pristine  vibration  data  is  available  anyway.  The  following  activity  indices  for  the  springy  with 
displacement  qa(t)  are  defined. 

ac  =7]ca(4  where  ca(t)  <  0  and  (  >  m  (7a) 

a,  =7Jca(Q|,  where  ca (t)  >  0  and  t  >  m  (7b) 
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The  subscripts  V  and  Y  refer  to  the  compression  and  tension  activity  in  the  spring  j.  The 
absolute  value  of  the  quantity  ca(t)  is  denoted  by  |c0(/)j .  The  bar  over  the  coefficients  ca(t)  is  an 
indication  of  the  pristine  vibration  database  relation, 


?o(0 
Hb  (0 


ca{t) 

ch{t) 


(8) 


Similarly,  these  activity  indices  for  the  springy'  with  displacement  qb(t)  can  be  inferred  by, 

Pc  =  Xl?a(0|;  where  caC0<0and  m  >  (!  (9a) 

A=^|co(0|,  where  ca(t)  >  0  and  m  >  t  (9b) 

Given  the  vibration  database  q(Y)  for  the  damaged  structure,  the  coefficient  ca(t)  is  computed 
following  the  transformation  in  Eq.  (3).  Hence  the  activity  indices  without  bar  for  the  damaged 


structure  are  computed  as  follows. 

<XC  =  EM0|,  where  ca(t)  <  0  and  t  >  m  (10a) 

a,  =  X  k  W|,  where  ca  (r)  >  0  and  g  >  m  (1  Ob) 

Pc  =  X k (t%  where  c„(0<0and  m>t  (1 0c) 

P,  =J]k(0|,  where  ca(r)>0and  m>£  (lOd) 


The  net  tensile  and  compressive  displacement  activity  index  (AI)  at  springy  is  defined  as, 

A!  =  Xfc -a’ |  +  £|a/  "«/|  +  +  Z|^  'A']  (U) 

i  !  i  i 

Ihe  superscript  V  is  an  indication  of  the  vibration  database  due  to  various  initial  conditions. 


Note  that  for  pristine  structures,  ac  =  ac ,  a,  =  a, ,  J3C  =  J3C  and  A  =  A  ■  Thus  AI  =  0.  When 

the  damage  is  present,  however,  AI  &  0 .  This  feature  holds  even  at  springs  where  damage  is  not 
present.  At  the  damaged  spring,  however,  A!  will  be  maximal.  This  property  of  AI  is  used  to 
localize  the  damage  in  the  springs.  For  the  springs  connecting  the  fixed  ends,  AI  is  computed  in 
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similar  fashion.  For  instance,  consider  the  first  spring  in  Figure  1.  This  spring  takes  compression 
when  q\(t)  is  negative  and  tension  when  q\(t)  is  positive.  Thus  define, 


Yc  =  X  ki  (f >1’  where  ch  (0  <  0 

(12a) 

Yc  =  W|.  where  9l  (0  <  0 

(12b) 

7,  where  9i(0>0 

(12c) 

F,  =  2>(0|,  where  4i(t)>° 

(12d) 

The  AI  for  spring  1  will  be, 

AI  =  XK'-^  +3;7'  ~r‘\ 

(13) 

i  l 


Note  that  spring  4  takes  compressive  load  when  q$(f)  is  positive  and  tensile  load  when  q$(t)  is 
negative.  Accordingly,  the  AI  in  Eq.  (13)  is  computed  using  the  following  steps. 


7  c  =  Xl^(0|. 

where  q3  (t)  >  0 

(14a) 

fc  =  2>(4 

where  q3>(t)>0 

(14b) 

il 

M 

Vi 

where  <73{0<G 

(14c) 

^1 

II 

M 

OJ 

V 

where  Z/3(/)  <  0 

(14d) 

In  the  next  section,  the  identified  damaged  spring  is  considered  for  growth  retardation. 
Depending  upon  the  spring  assumed  for  damage  growth  retardation,  a  mathematical  model 
suitable  for  LQR  design  is  presented.  Note  that  LQR  needs  all  state  variables  (displacements  and 
velocities)  as  measurements.  An  alternative  approach  would  be  to  use  observers  [30,31]  where 
these  state  variables  are  estimated  depending  upon  the  sensor  outputs. 

3.  Damage  growth  retardation  using  LQR 

Consider  the  transformation  given  in  Eq.  (5)  and  (6)  for  a  generic  pristine  structure  as. 
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q(r)  =  Hc(0 


05) 


such  that 


q(0  =  Hc(/)  and  q(f)  =  Hc(/) 


(16) 


A  finite  element  model  suitable  for  LQR  design  in  first  order  from  is  presented  in  sequel.  Recall 
the  model  of  the  pristine  structure  compatible  with  the  spring  and/or  damper  damages  as, 

Mf(r)  +  Cq(0+  K0q(0  -  Du(f)  ,  q(0)  *  0  . 

At  time  t  -  0 ,  let  q{0)  be  the  nonzero  initial  condition.  Without  loss  of  generality,  nonzero 
velocity  vectors  can  also  be  taken  as  an  initial  condition.  u(r)  is  the  control  input  vector  with  m 
components  to  be  determined  for  damage  growth  retardation.  Substituting  Eqs.  (15)  and  (16)  in 
the  pristine  model, 


M  H  c(f)  +  C0H  c(t)  +  K0Hc(/>  =  Du{/) 


(17) 


Eq.  (1 7)  in  first  order  from  is  given  by. 


—  (MH)-iKqH  -(MH)'’C0H 


(MH)"1  D 


u(0 


(18) 


Eq.  (18)  is  the  standard  structural  model  in  first  order  form  where  the  tensile  and  compressive 
coordinates  are  considered  for  the  cost  function  of  the  LQR.  One  of  the  important  concerns  in 
this  context  is  to  select  the  actuator  placements.  For  transient  vibration  monitoring,  any  actuator 
placement  would  work  as  long  as  it  is  used  to  place  eigenvalues  that  minimize  the  damage 
growth  retardation.  However,  for  steady  state  vibration  monitoring,  actuator  placement  becomes 
an  important  problem  as  the  actuator  load  is  expected  to  nullify  the  effects  of  structural  load  on 
the  damaged  spring  by  reducing  the  amplitudes  (effect  of  closed  loop  pole  locations).  In  practice, 
it  is  observed  that  the  transients  vanish  early  but  steady  state  vibration  is  equivalent  to  fatigue 
loads  with  constant  and  variable  excitation  frequency  spectrums.  Here,  it  is  expected  to  have  the 
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eigenvalue  distribution  of  the  damaged  structure  farther  from  those  of  the  undamaged  structure 
in  the  open  left  half  plane  of  the  complex  plane  so  that  fatigue  load  in  controlled  structure 
generate  less  vibratory  amplitudes  than  the  fatigue  load  in  uncontrolled  structure.  But  what 
eigenvalue  locations  minimize  these  effects  are  unknown  and  these  locations  may  vary 
depending  upon  the  number  of  actuators.  In  transient  case,  LQR  becomes  a  powerful  tool  to 
determine  such  eigenvalue  locations. 

If  the  model  in  Eq.  (18)  is  controllable,  then  the  LQR  design  procedure  is  recalled  briefly.  Let 
Eq.  (1 8)  be  of  the  form, 

i(0  =  Ax(f)  +  Bu(0  (19) 

The  control  forces  at  m  locations  that  minimize  the  cost  function  J  are  given  by, 
h(/)  =  -R'  'B  Px(0  =  -G  x(f)  (20) 

where  P  is  a  positive  definite  matrix  and  is  the  solution  of  the  algebraic  Riccati  equation, 
A’P+PA-PBR"lUP  +  Q  =  0,  (21) 

The  positive  semi-definite  matrix  Q  and  the  positive  definite  matrix  R  are  user  defined.  They 
define  the  cost  function  J  to  be  minimized  as  below, 


oo 

J  =  J(x'(OQxO)  +  u'(/)Ruft))rf/  (22) 

o 

For  the  discrete  structure  shown  in  Figure  1,  the  positive  semi-definite  matrix  Q  that  minimizes 
tensile  and  compressive  activity  at  the  damaged  spring  2  is  selected  as  follows. 

Case  A : 


10  0  0  0  0 
0  0  0  0  0  0 
0  0  0  0  0  0 
0  0  0  0  0  0 
0  0  0  0  0  0 
0  0  0  0  0  0 


if  c,2(0  is  to  be  minimized. 


(23a) 
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Case  B\ 


I  0  0  0  0  0 
0  0  0  0  0  0 
0  0  0  0  0  0 
0  0  0  10  0 
0  0  0  0  0  0 
0  0  0  0  0  0 


if  c12(0  +  c12(2)  is  to  be  minimized. 


(23b) 


For  simulations,  R  is  taken  as, 


R  =0.01 


1 

0 


0 

1 


In  case  A,  only  spring  2  is  assumed  to  be  damaged.  Whereas,  in  Case  B,  both  spring  2  and 
damper  2  between  the  masses  m\  and  m2  are  assumed  to  be  damaged.  Further  Case  A  takes  into 
account  only  the  potential  energy.  Whereas,  Case  B,  takes  into  account  both  potential  and  kinetic 
energies.  This  way  the  LQR  assists  to  retard  the  damage  growth.  However,  it  is  important  to 
observe  that  the  LQR  design  is  pursued  using  pristine  structure  and  it  is  applied  to  the  damaged 
structure.  Although  the  cost  function  is  theoretically  not  minimized,  it  is  expected  that  the  cost 
function  for  the  controlled  damaged  structure  is  lesser  than  the  cost  function  for  the  uncontrolled 
pristine  structure  if  the  eigenvalue  distribution  of  the  controlled  damaged  structure  is  farther 
from  those  of  the  controlled  pristine  structure  in  the  left  half  plane.  That  is,  the  eigenvalues  of 
the  controlled  damaged  structure  are  farther  from  the  eigenvalue  of  the  controlled  pristine 
structure  which  is  nearest  to  the  imaginary  axis.  This  is  equivalent  to  the  robust  alpha-stability 
problem  under  real  parameter  variations  which  depict  the  life  of  the  controlled  or  uncontrolled 
structures.  In  the  next  section,  admissible  stiffness  and  damping  parameter  variations  that  would 
retard  damage  growth  with  active  control  forces  determined  using  LQR  are  presented. 

4.  Admissible  damages  retarding  the  grow  th 
Consider  the  damages  A k2  and  Ae2  introduced  in  the  stiffness  and  damping  matrices  as  below. 


13 


K  =  K0  -  Ak2 K 2  and 


C  =  C0-Ae2C2  (24) 

K2  and  C2  are  the  known  constant  distribution  matrices.  A  fundamental  problem  in  damage 
assessment  for  structural  prognosis  is  to  know  the  magnitudes  of  the  scalars  A k2  and  Ae2  for 
which  the  eigenvalues  with  damaged  spring  and  damper  are  more  damped  and  less  oscillatory 
compared  to  those  eigenvalues  with  undamaged  spring  and  damper.  This  theme  applies  to 
controlled  as  well  as  uncontrolled  structures.  For  the  controlled  structure,  the  activity  index  is 
expected  to  be  more  if  the  damaged  structure  with  controller  exhibits  eigenvalues  closer  to  the 
imaginary  axis  and  away  from  the  real  axis  compared  to  the  undamaged  structure  with  the  same 
controller.  This  is  the  robust  alpha-stability  problem  where  alpha  is  the  eigenvalue  of  the 
undamaged  uncontrolled  structure  that  is  closest  to  the  imaginary  axis.  Note  that  the  farther 
imaginary  parts  of  the  eigenvalues  from  the  real  axis  contribute  to  the  frequency  of  oscillations, 
which  also  contribute  to  damage  growth  and  these  effects  of  eigenvalues  are  currently  ignored. 
The  robust  alpha  stability  problem  can  be  posed  as  follows.  Determine  the  marginal  A k2  and  Ae2 
such  that  the  eigenvalues  of  the  damaged  structure  with  controller  exhibit  at  least  one  of  its 
eigenvalues  on  the  alpha  line  parallel  to  the  imaginary  axis  and  the  rest  of  the  eigenvalues  are 
alpha  stable.  Mathematically,  these  conditions  are  formulated  as  follows. 

The  damaged  structure  with  controller  is  written  as, 

x(0  =  ((A  -  B G )  +  q{ E,  ■+■  #2  E  2  ) x(f)  =  ( A0  +  <7, E,  +  q2 E 2 ) jc(r)  =  A(qx , q2 ) x(/)  (25) 

Here,  the  matrix  A0  is  asymptotically  stable  with  its  eigenvalues  in  the  open  left  half  plane  of  the 
complex  plane.  The  matrices  E,  and  E2  compatible  to  qx  =A k2  and  q2  =Ae2  are  given  by. 


Ei  = 


0 

(MH)-'K2H 


o' 

0 
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(26) 


e2  = 


0  0 
0  (MH)_,C2H 


Admissible  damages  retarding  the  growth  are  the  admissible  q,,i  =  1,2  for  which  the  matrix 
A  {quq2)  is  alpha  stable  (or  a-  stable).  That  is,  the  system  with  matrix  \{qx,q2)  is  a -stable 
where  a  is  the  real  part  of  the  eigenvalue  of  the  matrix  A  which  is  closest  to  the  imaginary  axis. 
Also  observe  that  the  fatigue  cycles  of  compressive  and  tensile  activity  at  the  damaged  spring  are 
present  at  matrices  A(q: ,  q2 )  and  at  A  in  the  form  the  parameters  of  the  hat  function  in  matrix  H  . 
In  this  framework,  robust  analysis  tools  discussed  in  [22]  are  applied  to  determine  the  uncertain 
parametric  values  q,  for  a  -stability.  The  life  of  the  structure  when  control  is  applied  will  be 
extended  as  long  as  the  damage  growth  with  parameters  <?,  satisfy  a  -  stability.  These  parametric 
values  assist  to  monitor  the  life  of  the  structure  which  can  be  estimated  using  an  EKF  [10]. 

5.  Simulations 

In  order  to  illustrate  the  underlying  principles  of  structural  damage  localization  and  growth 
retardation  in  springs,  the  discrete  structure  shown  in  Figure  1  is  considered.  First  damage 
detection  in  spring  2  or  3  is  performed.  In  a  laboratory  setting  [32],  for  a  single  degree  of 
freedom  system,  it  has  been  suggested  that  by  measuring  amplitudes  and  frequencies,  the 
displacement  is  measured.  To  introduce  stiffness  loss  in  the  damaged  spring,  it  may  be  possible 
to  warm  up  the  spring  mildly  while  the  displacement  measurements  are  recorded.  But  for 
theoretical  investigations,  the  following  data  is  assumed  for  simulation.  mf  =  \,m2  =  2  and 
m3  =  3 .  kl  =  10 ,  k2  =  10 ,  £3  =  20  and  k4  -  20 .  The  damping  data  are  e)  =0.5,  e2  =  1 ,  e3  =1 
and  e4  -  0.5 .  For  damage  detection,  a  perturbation  in  spring  2  is  introduced.  That  is,  £2  is 
changed  from  1 0  to  9  (stiffness  loss  =  I  unit) ,  then  a  vibration  database  for  each  initial  condition 
set  (<7i(0)=0.1,  rest  are  zero;  c/2(0)=0. 1 ,  rest  are  zero;  <fc(0)=0.1,  rest  are  zero)  were  recorded. 
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Table  1  for  spring  2  and  Table  2  for  spring  3  and  Table  3  for  springs  1  and  4  summarize  the 
activity  indices  stated  in  Eq.  (7),  Eq.  (9),  Eq.  (10),  Eq.  (12),  Eq.  (13)  and  Eq.  (14)  respectively. 
As  discussed  before,  the  activity  indices  ac,  a,,  fic,  /?, ,  yc  and  yt  for  the  damaged  structure 

are  different  from  the  pristine  structure  ac,  a, ,  j3c>  j3 ,,  yc  and  y, .  Thus  Tables  1,  2  and  3 
indicate  damage  present  in  the  structure.  This  can  be  also  inferred  from  the  errors  a'  -  a'c , 

a,'  -  a\ ,  f5'c  -  fi‘c ,  j3't  -  p‘ ,  f'  -y'c  and  y\  -  y\  in  Figures  3  and  4,  respectively.  The  net 
activity  index  in  Eq.  (11)  and  in  Eq.  (13)  indicates  that  AI  for  Table  1  is  ‘2.2653’  and  for  Table  2 
is  ‘0.8739’.  Similarly,  the  net  activity  index  in  Eq.  (13)  indicates  that  Al  for  Table  3  using 
springs  1  and  4  is  1 .8666  and  0.5326,  respectively.  The  AI  against  the  springs  is  shown  in  Figure 
5.  Since  the  AI  for  spring  2  is  maximal,  it  is  concluded  that  spring  2  is  damaged. 

In  the  following  discussion,  damage  growth  retardation  in  spring  2  is  performed  using  the 
LQR.  Damage  assessment  criterion  for  growth  retardation  presented  in  Section  4  is  also 
illustrated  in  the  subsequent  discussion.  Figure  6  depicts  various  scenarios  by  which  the  transient 
fatigue  cycles  can  influence  damage  growth.  For  instance,  the  fatigue  cycles  of  an  uncontrolled 
structure  (in  red)  suggest  damage  growth.  Here  the  stiffness  of  the  spring  2  takes  the  nominal 
values.  Hence  it  is  categorized  as  no  damage.  The  transient  fatigue  cycles  of  the  controlled 
structure  (in  blue)  suggest  damage  growth  retardation.  Once  again  the  simulation  is  carried  out 
for  a  nominal  spring  stiffness  value.  Hence  it  is  categorized  as  no  damage.  Whereas,  the  transient 
fatigue  cycles  of  the  controlled  structure  (in  black)  repeats  damage  growth  retardation  scenario; 
however,  an  admissible  damage  size  for  growth  retardation  is  assumed  in  the  simulations.  The 
admissible  damage  size  is  determined  based  on  a  -  stability  criterion  discussed  in  Section  4.  The 
real  part  of  the  eigenvalue  closest  to  the  imaginary  axis  when  control  is  not  applied  is  indeed  a 
and  it  is  given  by, 
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a  =  -0.0831 1222  3797488  . 


Figure  7  depicts  the  a- stability  criterion.  The  figure  considers  uncontrolled  structure.  For 
the  stem  values  (admissible  damage  sizes  (stiffness  and  damping  loss)  for  growth  retardation), 
the  eigenvalues  are  preserved  with  a  -stability.  That  is,  for  these  values,  the  amplitudes  of  the 
transient  fatigue  cycles  with  damage  in  the  spring  are  retained  smaller  than  that  of  the  amplitudes 
oi  transient  fatigue  cycles  when  there  is  no  damage.  In  uncontrolled  structure,  although  this 
phenomenon  extends  the  life  of  the  structure,  it  is  further  required  that  the  admissible  damage 
sizes  for  growth  retardation  keep  the  structure  non-oscillatory.  This  requirement  depicts  the 
restrictions  on  the  imaginary  parts  of  the  eigenvalues,  which  are  difficult  to  formulate.  Also,  it  is 
observed  that  the  admissible  damping  coefficient  loss  for  damage  growth  retardation  increases 
with  stiffness  loss,  lhat  is  when  stiffness  loss  is  zero,  the  admissible  damping  loss  is  also  zero, 
^hen  the  stiffness  loss  increases,  admissible  damping  loss  also  increases  reaching  the  nominal 
value  when  M2  =  -7 .  In  the  controlled  structure,  this  property  is  preserved  for  all  stiffness  loss 
with  Ak2  =0,-1, -2**-,  etc.  Hence  the  admissible  damage  sizes  for  damage  growth  retardation  are 
indeed  the  values  at  the  stems  with  peak  values  approaching  the  nominal  values  of  the  spring  and 
damper.  Note  that  for  any  admissible  damage  size  that  might  take  in  real  time,  the  amplitudes  of 
the  transient  fatigue  cycles  are  always  smaller  than  those  of  the  undamaged  structure.  This 
phenomenon  also  applies  to  the  rate  at  which  transient  fatigue  cycles  change  with  respect  to  time 
as  shown  in  Figure  8. 

Next,  the  above  eigenvalue  based  damage  growth  retardation  discussion  is  validated  with  the 
time  responses  of  C|  (t)  and  c,  (?) .  In  Figure  9,  the  compression  activity  is  denoted  by  negative 
values,  these  values  are  sum  of  cY (/)  over  t  when  c, (7) < 0 .  Similarly,  the  tensile  activity  is 
denoted  by  positive  values.  These  values  are  sum  of  c,(r)  over  t  when  c,  (t)>  0 .  The  uncontrolled 
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structure  exhibit  Larger  values  compared  to  the  controlled  structure.  In  both  the  cases,  clearly, 
damage  growth  retardation  is  observed.  That  is,  when  stiffness  loss  increases,  the  displacement 
activity  also  increases.  Note  that  for  a -stability,  the  stiffness  loss  is  restricted  to  5.  In  these 
simulations,  the  damping  loss  is  assumed  zero.  Similar  observations  are  recorded  with  tensile 
and  compressive  activity  rates  in  Figure  10.  These  simulations  presented  above  suggest  the 
validation  of  the  technique  for  in  situ  damage  localization  and  damage  growth  retardation 
discussed  in  the  paper. 

6.  Conclusions 

Given  an  operational  structure  or  a  mechanical  system  modeled  as  a  spring-mass-damper 
system,  this  paper  presents  in  situ  damage  localization  and  assessment  with  sizes  and  location  of 
the  damage  determined  completely  for  damage  growth  retardation.  The  new  technique  for 
damage  localization  using  isolated  hat  functions  applied  to  a  vibration  database  is  presented. 
After  the  damage  is  localized,  damage  growth  retardation  at  the  spring  is  pursued  using  the  linear 
quadratic  regulator  design.  The  weighting  matrices  are  particularly  attractive  to  minimize  fatigue 
cycles.  But,  the  design  is  performed  by  using  a  pristine  structural  model.  Hence  it  becomes 
necessary  to  configure  damage  growth  retardation  when  the  model  is  damaged.  This  problem  is 
posed  as  a  robust  alpha  stability  problem,  where  the  damage  assessment  is  carried  out  to 
guarantee  damage  growth  retardation.  A  spring  mass  damper  system  is  considered  to  illustrate 
the  damage  localization  and  growth  retardation  problem  considered  in  this  paper. 
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Table  1  Load  Test  (Pristine/ (Damaged))  for  Passive  Damage  Detection  in  Spring  2 


Displacement 

<h  (0 

?2(0 

Activitylndex 

a,  /(a,) 

a,/(ac) 

A /(A) 

A,  /(A) 

£k„(OI 

3.9343 

4.6297 

4.0989 

4.5737 

9,(0)  =  0.1 

(3,9502) 

(4.6894) 

(4.0074) 

(4.5092) 

Zkwi  i 

5.7751 

5.2936 

12.2315 

11.7283 

o 

II 

o 

<N 

(5.3742) 

(4.8758) 

(12.4116) 

(11.8263) 

2>„<oi 

6.8082 

6.6592 

12.6969 

13.1245 

9,(0)  =  0.1 

(6.4400) 

(6.3169) 

(12.8291) 

(13.2187) 

Table  2  Load  Test  (Pristine/ (Damaged))  for  Passive  Damage  Detection  in  Spring  3 


Displacement 

(0 

?3(0 

Activity  Index 

«,  /(or, ) 

acl(ac) 

A  /(A) 

A /(A) 

2>„(0I 

4.0790 

3.8611 

2.7193 

2.7946 

9,(0)  =  0.1 

(3.9211) 

(3.6904) 

(2.6008) 

(2.6749) 

Zk.wi 

11.4170 

11.8512 

6.8322 

7.1069 

<7,(0}  =  0.1 

(11.4150) 

(11.9163) 

(6.8313) 

(7.0778) 

2>„(')l 

13.6455 

13.0964 

9.0454 

8.8299 

9,(0)  =  0.1 

(13.6116) 

(13.1019) 

(8.9710) 

(8.7337) 
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Table  3  Load  Test  (Pristine/(Damaged))  for  Passive  Damage  Detection  in  Springs  1  and  4 


Displacement 

<i\  (0 

tf3(0 

Activity  Index 

y,Kr,) 

rJ(ye) 

7,  Ky,) 

7c  !(Jc) 

2>,(oi 

6.0424 

5.4054 

4.1319 

3.9564 

i  i 

o 

II 

o 

S3* 

(6.0320) 

(5.3537) 

(3.9857) 

(3.8066) 

Zktoi 

9.8675 

10.2396 

1 1 .6773 

11.5232 

<?2  (0)  =  0.1 

(9.4312) 

(9.7980) 

(11.6182) 

(11.5240) 

Zk(')i 

12.0426 

12.3883 

13.9516 

13.9901 

?3(0)  =  0.1 

(11.5874) 

(11.8873) 

(13.8425) 

(13.9225) 
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Errors  Errors 


3  4  5  6  7  8  9  tO  11  12 

Cell  Numbered  Row  Wise  {Table  1) 


0.3  | - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - r 


_q  - 1 _ I _ i _ i _ i _ i _ i _ i _ i _ i _ i _ I _ 1— 

1  2  3  4  5  6  7  8  9  10  11  12 

Cell  Numbered  RowWise  {Table  2) 


Figure  3  Errors  Contributing  to  Activity  Index  in  Eq.  (1 1) 
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Errors 


Ceil  Numbered  Column  Wise  (Table  3) 

Figure  4  Errors  Contributing  to  Activity  Index  in  Eq.  (13). 
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Figure  5  Damage  Localization  Using  Maximal  Activity  Index  (AI)  at  Spring  2  lor  a  Given 
Vibration  Database  From  the  Structure  in  Figure  1 . 
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28 


Admissible  damping  loss 
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Figure  8.  Uncontrolled  (red)  and  controlled  (blue  and  black)  fatigue  rates  favoring  (red)  and  retarding  (blue  and 

black)  the  damage  growth  in  spring  2, 
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Tensile  and  compressive  activity 


D _ i _ _ I - - - 1 - J - 1 - L - 

'6  0  1  2  3  4  5 

Stiffness  loss 

Figure  9:  Robust  alpha-stability  criterion  for  damage  growth  retardation  validated  by  displacement  activity. 
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Tensile  and  compressive  rates 


Stiffness  loss 

Figure  10:  Robust  alpha-stability  criterion  for  damage  growth  retardation  validated  by  velocity  activity. 
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